Moil. Not. R. Astron. Soc. OOP. HUMl (2002) Printed 12 January 2013 (MN 1^1^ style file v2.2) 



Profiles of spectral lines from failed and decelerated 
winds from neutron stars and black holes. 

A. V. Dorodnitsyn^* 

^Laboratory for High Energy Astrophysics, NASA Goddard Space Flight Center, Code 662, Greenbelt, MD, 20771, USA 



ABSTRACT 

We calculate profiles of spectral lines from an extended outflow from the 
compact object (a black hole, or a neutron star). We assume that the bulk 
velocity of the flow increases during a short phase of acceleration and then 
rapidly decreases forming a failed wind. We also study the wind which is only 
decelerating. We show that depending on the relative strength of the gravita- 
tional redshifting, line profiles from such winds may be of several types: dis- 
torted P-Cygni (emission and blueshifted absorption); W-shaped (absorption- 
emission-absorption); and inverted P-Cygni (emission - redshifted absorption). 
The latter case is expected from accretion flows where the velocity is directed 
inward, however we show that inverted P-Cygni profile can be produced by the 
failed wind, provided the line is formed within several tens of Schwarzschild 
radii from the compact object. 

Key words: line formation - radiative transfer - galaxies: active - radiation 
mechanisms: general - stars: mass loss - stars: winds, outfiows 



1 INTRODUCTION 



Relativistically broadened, double horn profiles from fluorescent Fe Ka lines is a widely 
accepted evidence for the presence of the cold accretion disk in Seyfert 1 active galactic 
nucleus (AGN), and the assumpti o n tha t such disk is illuminated by X-rays generated by 
the corona (see, e.g. iFabian et al.l ( l2000l )). However, the geometrically broad line forming 



region spans a large range of orbital velocities which translates to a broad region in frequency 
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space, additionally skewing the line by relativistic effects and making detailed diagnostics 
of the accreting plasma difficult. 

Spectral resolution of grating spectrographs of the X-ray telescopes Chandra and XMM- 
Newton below ~ 10 keV all ows observations re vealing complicated structure of the Ka line, 
in m any Type 1 AGNs e. g. (IReeves et al.ll200ll ): and for the detection of the narrow core, see 
Yaqoob et al.l ( 120031 ). Observations suggest that distant parts of the accretion disk as 



also 



well as the reflection from the obs curing torus are involved. Multi dimensional simulations of 



spectra of warm absorber flows by 



Dorodnitsyn fc KallmanI (120091 ) also show the importance 



of distant AGN winds for the formation of the complex structure of Fe Ka line. 

It is widely accepted that strong gravitational field of the compact object, such as black 
hole (BH), or neutron star (NS) can imprint itself into the radiation which comes from 
regions within several tens of Schwarzschild radii, Vg. Narrow emission and absorption lines 
from black holes and neutron stars delivers an optimal opportunity to deduce their masses 
and radii (in case of NS). 

Narrow lines formed near a BH carry direct information about the dynamics, ionization 
stage and covering fraction of the absorbing gas and thus about the mass budget of the 
inflow/outflow. The latter is important for understanding of the fraction of the AGN - host 
galaxy feedback which comes directly from regions close the BH. Gravitational redshifting 
is entangled with Doppler shifts, and altogether they encode information about the mass of 
the BH. 

Observed narrow absorption features from several Quasars and Seyfert galaxies possibly 
reveal the importance of gravitational redshifting and Doppler blue-, and redshifting due 
to the bulk motion of the moving gas. Such evidence was found in several cases includ- 
ing: Seyfert 1 galaxy NGC 3516 where it is suggested that observed absorption features 
ar e from the gravit ationally redshifted resonant line scattering within the accretion flow 



(INandra et al. 



19991 ): the gravitationally redshifted (10 - 20 Vg from the BH) resonance ab- 



sorption line from FeXXV or FeXXVI in the quasar E182 1+643 (lYaqoob et al.ll2005l ): a 
combination of the gravitational and Doppler shifting of lines fror n ionized Iron i s sugg ested 



to explain absorption features from the quasar PG 1211+143 by 



Reeves et al. 



see 



Matt et al. 



(120051); also 



( 120051 ) for possible detection of the strongly redshifted, transient absorption 
fea ture from ioni z ed Fe in the quasar Q0056-363. 

Dadina et al.l (120051 ) reported on the redshifted absorption, and emission features ob- 



served at ~ 6 keV, from Seyfert 1 galaxy Mrk 509. Two of the absorption features are 
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separated by the emission line, and the overall emission structure is attributed to H-, or 
He-like Iron. This W-shaped structure is found to be sporadic with variability as short as 
~ 20 ks. Infall and outflow, ~ ± (0.1 — 0.2) c of matter are suggested, together with the 
gravitational redshifting to explain these features. Transient nature of such features is con- 



sistent with the latter non-detection of t 



blueshifted absorption lines ( ICappi et al. 



le red shifted absorption and the detection of only 



20091). 



The detection of the redshifted absorption lines form accreting neutron stars which ex- 



putting 
(I2OO2I) 



hibit thermonuclear X-ray bursts are ideally suited for deducing the NS mass and 
constrains on its equation of state. So far, these detection are rare. ICottam et al. 
identified several absorption features in the burst spectra of the neutron star EXO0748-676. 
Lines from H-like Iron for the early phases of the burst, and He-like Iron for the late phases 
were adopted to explain the spectra, and to obtain the redshift oi Z = 0.35. Other stud- 
ies, including LTE and non-LTE neutron star atm osphere mo deling suggested n = 2 — 3 



transition of FeXXIV, and the redshift Z = 0.24 (iRauch et al 



20081 ). These lines are ex- 



pected to be intrinsically transient which additionally lowers the statistical significance of 
such detections. 

Usually, the existing models of gravitationally redshifted lines do not include influence 
of the transfer effects in the extended moving envelope, and if they do, simple transmission 
models are adopted. Large intensities of emission components, and simultaneous presence of 
emission and absorption from the same line transition suggest the importance of the coupling 



between dynamics of the win d and g ravitationa 
Beginning from works of 



redshifting. 



Bealsl (11929. 193ll ) it was understood that P-Cygni profiles 



provide a unique evidence for the rapidly moving wind. It was also understood that radiation 
from a huge volume, occupied by such wind naturally explains bright emission features in 
the spectra, and the observed broadening is due to high velocity dispersion of such wind 
(few X lO^kms"^ in the case of a normal star). The most remarkable about this profile is, 
of course, the relative blueshifting of the absorption trough relative to the emission line. As 
a P-Cygni profile bears an imprint of the bulk motion of the plasma, it is an important tool 



of the diagnostics o 
In the paper by 



the dynamics of st ellar winds. 



Dorodnitsyru (l2009l ) (hereafter Paper I) a modified Sobolev approxima- 



tion was built which incorporate gravitational redshifting, and spectral line profiles were 
calculated from a stellar-type and explosion-type winds, accelerated in the strong gravita- 
tional field of the compact object. 
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An important parameter from Paper I measures the relative importance of gravity at the 
wind base go = Rc/vg, where -Rc is the radius of the wind launching point. Two types of the 
velocity laws were adopted: i) "stellar" type, i.e. v{r) ~ (1 — Rc/r)^, and ii) Homologous 
expansion i>(r) ~ r (Hubble law), adopted for the description of outbursts. Important also 
found the terminal velocity, V°°, which spans a range of 0.01- 0.3 c and the distribution of 
the opacity. 

It was found that gradually accelerated winds which are launched at i?c ^ 20 — 30 
produce the following types of profiles: i) distorted P-Cygni, i.e. characterized by emission 
feature which is redshifted relative to an absorption trough ii) saw-tooth, i.e. the narrow red- 
shifted absorption line superimposed on a broader redshifted emission line, combined addi- 
tionally with a blueshifted absorption line iii) W-shaped, i.e. absorption-emission-absorption 
profiles (by "red-shifted" or 'blue-shifted" we necessarily mean only the relative blue-, red- 
shift of one spectral feature with respect to the other, which may, or may not represent the 
actual shift of frequency of such feature with respect to the rest frequency of the line). 

In this paper we calculate profiles of lines from a wind which is either everywhere decel- 
erated, i.e. described by the i> ~ law or failed, i.e. its velocity rises steeply during the 
short acceleration phase, and then decreases, asymptotically approaching v ~ law. 

We will show that fine profiles from failed winds are significantly different from those 
considered in Paper I. Additionally to three types described above there is a forth type: the 
inverted P-Cygni profile (absorption-emission). Notice, that the inverted P-Cygni profile 
is expected from spherically-symmetric accretion fiow, and to observe it from the wind 
is counter-intuitive. We will show that strong gravitational redshifting together with the 
"failed" character of the wind, and also provided the deceleration phase is fast (i.e. m > 3) 
are necessary conditions for the appearance of this feature. Failed wind profiles can also be 
W-shaped or saw-tooth shaped. 

The plan of this paper is as follows: in Section 2 we review basic assumptions of Paper 
I about the wind geometry and gravitational field; in Section 3 we describe our methods of 
calculation of fine profiles i.e. the Sobolev optical depth of the fine, source functions, calculate 
mean radiation field etc.; in Section 4 we describe velocity laws, derive the bridging formula 
for the failed wind and make assumptions about the distribution of the opacity; in Section 

5 we calculate equal frequency surfaces, which arc of paramount importance to understand 
line profiles which we calculate in Section 6; in Section 7 and 8 we discuss major results and 
conclude. 



P-Cygni profiles from outflows near compact objects. 5 
2 ASSUMPTIONS AND APPROXIMATIONS 

We assume that photons from an external source of continuum radiation are interacting 
with the wind via resonance hne scattering; both the wind and the radiation source are 
assumed to be spherically symmetric; the wind is launched from the photosphere of the 
radius Rc\ there is a prescribed velocity profile, v{r) and a prescribed distribution of the 
line opacity in the wind. The only difference in the assumptions with those of Paper I is 
that here we consider decelerating winds, instead of winds which are gradually increasing 



their velocity. Profiles 
different authors (e.g. 



Tom decelerated winds from normal stars were ca^ 



Kuan fc Kuhil ( 119751 ): Marti fc Noerdlingerl (119771 )). In the following 



culated in works of 



we briefly summarize the formalism developed in Paper I. 

We assume that after leaving the photosphere a photon is traveling in the wind without 
interaction with the matter except for certain points (resonances) where it's frequency in 
the co-moving frame, z/ appears to be within the Doppler thermal width of the spectral 
line, Az//) = v^vth/c, where Vth is the thermal velocity, and z/q is the frequency at the line 
center. To allow for the resonance region to be located as close as several Schwarzschild radii, 
Tg = 2GM/c^ from the compact object we take into account both Doppler and gravitational 
shifting of the photon's frequency. The Lorentz transformations from a local Lorentzian 
frame (a photon frequency z^ioc) which is at rest at a given point Sq to the co-moving frame 
give: u = 7z/ioc(l — /i/3), where /3 = f/c, 7 = (1 — /3^)~^/^. A photon that was emitted at the 
point So after traveling to some other point s is gravitationally red-, blueshifted and in the co- 



moving frame has the following frequency: z/(s) = 7z^ioc('5o) V^oo(^o)7^oo(^ (1 ~ 

Gravitational field is described in a "weak field limit", i.e. by means of the gravitational 
potential, 0, and also implying y^goo — 1 + To allow for the terminal velocity, to 
be of the order of the escape velocity V^sc at the base of the wind (which can be a large 
fraction of the speed of light, c), we must retain all terms of the order of v'^/c^ and 

The photon emitted at the point s, in a direction n, having the frequency z> will be 
registered by the observer at infinity at the frequency z/°°: 

sW = .~(i-,W/JW-^ + ^). (1) 

where /z = n ■ v = cos^^, and is the gravitational potential. 

It is convenient to adopt dimensionless units: the non-dimensional frequency, y = {u — 
z/o)/Az/D, the non-dimensional radius, x = r/Rc, the non-dimensional velocity, u = v/vth, 
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and the non-dimensional gravitational potential $ = (j)/{cvi\^) = C(0/c^), where ( = c/vth- 
In the narrow line limit, the emitted frequency is uq and thus 

= (zy- _ ^q)/Az/d. (2) 

In such non-dimensional units, the equation ([T]) can be cast in the form: 

y^=fiu{x) + ^x)-u{xY/{2C). (3) 

Equation determines the position of the resonant point. It is straightforward to rewrite 
the above equations in terms of the redshift, Z = {u — z/°°)/z/°°: 

y- = -Cjf^, (4) 



and from equation ([3]), we obtain: 

/iM + $-nV(2C) 



(5) 



C + /iM + $-«V(2C)' 

To describe gravitational field we adopt two types of potentials: the Newt onian potential, 
(f)(r) = —^^^ and the pseudo-Newtonian potential of Paczynski-Wiita (PW) (jPaczynski fc Wiita 



20051): 

GM 



Tg — r 



(6) 



The latter mimics important features of exact general relativistic solutions for particle tra- 
jectories near a Schwarzschild black hole. For example, it reproduces the positions of both 
the last stable circular orbit, located at 3rg and the marginally bound circular orbit at 
2rg. The non-dimensional form of the PW potential, ([6]) reads: = (/ (2(1 — xgo)). A 

non-dimensional parameter, qq = Rc/^g determines the relative importance of the gravi- 
tational redshifting (i.e. by equating ^ oo one completely neglects the influence of the 
gravitational field on the energy of a photon). 



We adopt a {p, z) coordinate system, such that x = \/ p"^ + z"^, where p is the impact 
parameter and the observer is located at 2; = oo. In the non-dimensional variables, equation 
([3]) takes the form: 

zq I r~o, \ 9 \ , I ( r~o, \ o \ 1 



The solution of the equation ((Tj) for different values of p determines the locus and shape 
of the equal frequency surface (hereafter EPS). Understanding the shape, and topology of 
EPS is very helpful in order to understand and interpret profiles of spectral lines from the 
moving medium. Notice, that equation, ([7]) does not reference to any particular form of the 
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velocity law, v{r). In the absence of gravitational redsliifting, the first term on the right hand 
side of (I7j) determines the locus of the surfaces of equal line-of-sight velocities. The second 
and third terms, stand for the gravitational redshifting and transverse Doppler effect (also 
redshifting) . Equation ([7]) was derived in Paper I, where its properties are investigated for 
the case of the outwardly accelerated wind. Here we study this equation from the perspective 
of decelerated winds. 



3 CALCULATION OF LINE PROFILES: METHODS 



To include gravitational redshifting into the frame of the Sobolev approximation, requires 
to take into account terms of the order (9(f^/c^). Sobolev optical depth enters Castor's 
( ICastorl Il970l ) escape probabilities, thus not only making them dependent on d(f)/ds, but 
also demanding to i nclude 0(v^/c'^) terms int o their derivation. A Sobolev approximation 



e.g. (jSobolev 



1960; 



Rybicki &: Hummer 



1978f) which takes into account all of the above 



assumptions, was developed in Paper I, and including decelerating wind into consideration 
leaves this formalism intact. The Sobolev length, Lgob is defined as a characteristic length 
over which the photon, propagating in the direction s stays in resonance with a certain line 
transition: Lsob — Vtt/ \dv/ds + (l/c) d(f)/ds\. The optical thickness in the direction s is 



X L 



sob- 



Tsoh — Xo.com 

The Sobolev optical depth in a line is found from the relation: 



(8) 



n 



Xo ,com 
1 



rfln/3 
dlnr 



+ /i 



/3 



d\n(3 



/3 c'^ dlnr dlnr 



where Xo com opacity at the line frequency in the co-moving frame: 



Xo, 



vre 



(gf) 



Ni/gi - Nu/gu 



(9) 



(10) 



mc '~ " ' Az/ 

where A^^, Ni and gi are populations and statistical weights of the corresponding levels 
of the transition, and / is the oscillator strength of the transition. Retaining terms of the 
order C(t>^/c^), we recast equation ^ in the form: 



Tr X 



dlnB , ^ , , , 1 

^l + /^(l-2/i)) + 



dlnr 



c^f3 dlnr 



(11) 



8 A. V. Dorodnitsyn 



X 



1 - 



d\n.r 



-(3 



d\nl3 



1 -1 



JS d\nr d\nr 

where Tr = ri(yU = 1) is the optical depth in the radial direction. The Sobolev approximation 
consists of two parts: i) calculation of the optical depth of a resonant layer (i.e. from (j9])), and 
ii) calculation of the source function 5*1,. Considering the scattering within a line and adopting 
complete frequency redi stribution , and i sotropic scattering matrix, the source function can 
be approximated as (i.e. iHummerl (jl969[ )): 



(12) 



where Ji^ is the mean intensity, and Bi, is the Planck function, e is the ratio of coUisional 
and de-exitation rates. 

In order to calculate the source fun ction and the n, finally, the emergent intensity, we 
adopt a method of escape probabilities of lCaston (119701 ). In this approach, the mean intensity 
Ji/ is found from averaging of the formal solution of the radiation transfer equation over all 
solid angles; the result is 

Ju = Su{l — Pesc) + ^dist(-Ppen); (13) 

where Pesc is the probability for a photon to escape from the resonant region, and Jdist 
is the contribution to the mean radiation field delivered by the distant core, and Ppen is 
the probability for a photon, which is issued by the core, to penetrate to the given point. 
Combining (fT2l) and (|T3l) in a single linear equation and solving for 5*1^, we obtain: ~ 
<^ist (-Ppen) /-Pesc (uoticc, that wc assumcd e = 0). The probability for a photon to escape in 
a direction /i reads: 

Pesc = {PL (1 + /3'(3/i' - 1) + 2/i/3)), (14) 

where () denotes angular averaging: {f)n = ^ f /(^') dQ', and Pl^^ is the non-relativistic 
expression for the directional escape probability: 



-'loc 



1 - exp(-ri(/x,s)) 



(15) 



ri(s,/i) 

where ri is obtained from (Q. We assume that the core emits continuum black-body radi- 
ation, with no limb darkening, and that in the vicinity of the line the continuum may be 
approximated as constant. Integrating Pl^^ over the maximum angle subtended by the core: 



9c = arccos /Xc = arccos a/1 — x^^, the probability for a photon to penetrate to a given point 
reads: 



pen 



(PL) 



(16) 
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where is a dilution factor: 



The distant contribution, J^ist is calculated from 



•-^dist -^c Ppeni 



and for the source function, Su we have 
J. 



(l + /32(3/i2-l) + 2/i/3) 



(17) 



(18) 



(19) 



p 

esc 

In the presence of non-monotonic velocity law, for example when the wind first accelerates 
and then decelerates, multiple resonances may occur. This are best studied adopting a 
concept of equal frequency surfaces (see next section). At the i^^ resonance, the radiation 
field suffers a change: some photons from directions other than that of the ray, are scattered 
into the direction of the observer (local contribution); the other contribution is due to the 
incident intensity which is attenuated at this point: 



(20) 



where Si is the source function at the i*^ resonance. The source function Si is anisotropic 
in the direction of m o tion a nd was calculated in the relativistic case in the paper by 
Hutsemekers fc Surdejl (1l995l ). Notice, that these authors calculated the full relativistic 
source function, and our calculations take into account the contribution from the gravita- 
tion shifting, and also from the relativistic effects which are included by retaining (9(f^/c^) 
terms. 

We assume that resonant interaction of photons and matter takes place while propagating 
only in the forward direction (in the direction towards the observer), and that all other pos- 
sible interactio ns are neglected. Altogether, these assumptions are known as the disconnected 



appro ximation (j Grachev fc Grinin 



1975 



Marti fc Noerdlinger 



1977 



Rybicki fc Hummer 



19781 . 



19831 ). It is usually adopted in the multiple resonance case (see Paper I for the discussion of 
approximations) . 

Our computational domain has a spherical shape of a radius Pout- A ray with impact 
parameter, p intersects this sphere in two places: zi, and Z2 > zi. To obtain the emergent 
intensity, we sum the expression ( 120|) over all resonances along the ray with impact parameter 
p. The total number of resonances is not known a priori and thus the procedure of the 
calculation of the total intensity is best cast in the form of a simple recurrence formula: 
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h{v°°,p) = /._ie-^'- + 5,(1 - e-^>-), (21) 
where i = 1..N, and is the total number of resonances, the initial intensity Ii reads: 

where the subscript refers to the point which is the farthest point along the ray from the 
observer in the computational domain i.e. closest to Zi. Notice, that f l2ip is the summation 
formulae, it accounts for the fact that we don't know the total number of resonances before 
we transfer a photon through the computational domain. Thus we propagate the photon 
along it's trajectory (i.e. along a straight line with impact parameter p, and checking for 
the resonance with the line. At the i*'* resonance we calculate Jj , and go further; if another 
resonance is detected, this, previous Jj becomes a in a new iteration. Integration along 
the photon's trajectory is performed until the point Z2 is reached. All other couplings between 
resonances are neglected. Non-local coupling b etween resonance s in the case of a normal hot- 



star, non-monotonic wind was investigated in lPuls et al.l (119931 ). From their results we may 



expect our approximation of the photon transfer is accurate within 10%. 

After being normalized to the continuum flux, Fc the radiation flux seen by the observer 
is calculated from the following: 



00 



F(z/-,p)/F,= / I^{u^,p)pdp, (23) 
Jo 

where = /Ar(z/°°,p). 



4 DISTRIBUTIONS OF THE VELOCITY AND OPACITY 

In this paper we consider two types of the velocity laws. Both are used to approximate a 
wind which failed to escape from the potential well. The flrst adopted law describes the 
decelerating wind, 

u{x) = U-^, (24) 

where m > 0, and U°° = V°°/vth, V'^ is a terminal velocity. Equation describes the 
wind, decelerating from r = where its velocity has a maximum. A different approxima- 
tion is used to describe the wind which is accelerating at the beginning but after reaching 
maximum velocity is decelerating: 

uix) = Uo(^-il- wo/Uo)e-^^-'^"'A (25) 
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where Wq = v{Rc)/V°° = 0.01 and < e < 0.1 is related to the location or "thickness" of 
the boundary layer, which allows a smooth matching of the inner, increasing solution with 
the outer, decreasing one. The value of Uq is related to U°°, and is calculated numerically so 
that ma.x{u{x)) = U°°. In the acceleration part, in the limit 5x = x — 1 ~ 0, equation (!25|) 
reduces to u{x) ~ const ■ {6x)"^ + Wq, and at infinity, l/x"^ term in the equation fl25|) is the 
dominant one. Velocity profiles and (125|) are shown in Figure [1] (left panel) for different 
values of (m, e). In the Paper I, explosion- type outbursts (described by the homologous, 
Hubble- type velocity law), and winds with a stellar-type velocity profile were considered. 
With simplifications already adopted in this paper one can argue that there is no need in any 
further specification of the particular physical mechanism which could have lead to velocity 
laws i^^, or (!25l) (although see discussion in Paper I). Having said that we consider them 
as a mere parameterization which may be releva nt to the real failed or de celerated winds. 

The distribution of the opacity is adopted from lCastor fc Lamerd (119791 ). who made use of 
the fact that in the optically thin case and the photoionization by diluted stellar continuum, 
the ionization ratio of the successive stages is proportional to W/n, where W is the dilution 
factor and n is the number density. With the help of the continuity equation, in a spherically 
symmetric geometry, this ratio is proportional to f . It results in the parameterization of the 
radial optical depth, Tr = ri(yU = 1) in terms of the non-dimensional ve locity, w = v/V°°, 
and is written in the following form: Tr ~ w. , Castor fc Lamerd (119791 ) also argued that 
Tr ~ 1 — 'u; law is somewhat more appropriate as the resultant P-Cygni profiles are closer 
to the observed ones. Notice that in the absence of self-consistent modeling of the opacity 
distribution together with the dynamics of the flow any such profiling of the opacity should 
be considered with caution. Most important for our means is weather Tj. peaks close to the 
photosphere or at large radii.. We chose the following parameterization of Tj.: 

r,(w)=To'^{l-w'), (26) 

where the parameter, k ^ 0, and the parameter Tq is related to the total optical depth at 

the line center, / dw = Tq. Figure [T] (right panel) shows distributions of the opacity (12^ 
Jo 

calculated from velocity laws (|24|) . ( 125|) for different pairs of (m, e), and for k = 1 and k = 2. 



5 SURFACES OF EQUAL FREQUENCIES 

Equation ([3]) determines the locus of equal frequency surfaces (EFS). That is, for a particular 
choice of the impact parameter, p and the frequency, y°° this equation determines the position 
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Figure 1. Velocity (left) and opacity (right) laws. Velocity curves: thin solid line: equation l|24ll : thick solid line: equation 
I I25II : curves are labeled by pairs of m, e. Opacity curves calculated from the equation II26I I for corresponding velocity profiles; 
thin solid line: using equation I I24I1 . and thick solid line: using equation II25I I: curves are marked by pairs of parameters p, and 
e from equation II26I I 



of the resonant point, Zq. In the Paper I it was shown, that in strong gravitational field, EPS 
have complicated shape and topology. Those of them which are calculated for the redshifted 
photons, may have branches situated between the star and observer. Of course, in the case 
of a normal stellar wind photons which are emitted in the gas which is approaching the 
observer, are blueshifted in the observer frame. 

In the language of resonant surfaces, it means there are no redshifted EPS in front of 
the star (we call redshifted or blueshifted EPS(z/°^) depending on whether is blue-, or 
redshifted). If gravitational redshifting overwhelms Doppler boosting then in the observer 
frame these photons are redshifted. 

It is useful to imagine resonant surfaces as partially transparent screens, which can scatter 
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photons of a certain frequency. If EFS situates in front of the star, it produces an absorption 
hne, according to /oe~^', where Jq is the incident intensity and ti is the hne optical depth, 
otherwise, EFS is reflecting, i.e.: / = (1 — e~'^') S, where 5* is the source function. Notice 
that EFS also can be multi-branched. 

The topology and shape of resonant surfaces in the case of the decelerated wind are 
different from those found in Paper I for the case of monotonically accelerated wind. In 
Paper I it was found that sufficiently strong gravitational redshift can overwhelm doppler 
blueshift which results in a redshifted EFS positioned in front of the star. 

To calculate resonant surfaces we introduce the following non-dimensional parameters: 
A = ?/°°/?/max — y°°/uraa.x, whcrc Mmax = U°° in the casc of the velocity law, 1^^, or the 
maximum velocity calculated from the velocity law, fl2^ . and the additional parameter 
/^max = (^imax/c)fth- Rewriting equation ([7]), we obtain: 

A = fiw{x) + — ^ - w{xfP^,,^. (27) 

2(1 - Xgo) /3max 

Equation ( 1271) is equivalent to equation ([T]). Notice that if no relativistic corrections are taken 
into account (i.e., no last term on the right), and additionally, the Newtonian potential is 
assumed, the second term on the right becomes — l/(2(yfo /3max) = const x x~^, and the 
problem, is described in terms of only two parameters. 

Equal frequency surfaces calculated for the velocity law flM|) . are shown in Figure [2l 
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Increasing gravitational redshift, significantly changes the shape and position of EFSs, 
advancing the redshifted EFS in front of the star, and leading to formation of the redshifted 
absorption line or absorption edge. 

Figure [3] shows EFS for velocity laws ( 12^ and ( l25l) for go = 16, m = 2. There are 
branches of both redshifted, and blueshifted EFS in front of the star, which makes possible 
the formation of both red- and blueshifted absorption lines (W-shape profiles). From Figure 
[3l we see, that EFSs of the failed wind are not very different from those of the decelerated 
wind. 

Each node in the {pi} grid marks a ray, and the problem of calculating the line profiles 
in p, z geometry is reduced to summation of contributions from all such possible rays. Each 
such ray is discretized, {z^}, and we iterate from zi to Z2 looking for resonances, and solving 
if necessary the non-linear equation ([7]). 
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Figure 2. Equal frequency surfaces for the velocity law l|24| l. The parameters are: go = 2000, /3max = 0.01 {upper left), 
go = 10, 13 max — 0.1 {upper righ£), go — 8.33, jSrYiax = 0.05 {lower left), go = 5, P max — 0.02 {lower right). Curves are marked 
by the parameter A, corresponding to the non-dimensional frequency at infinity. Curves: thick solid line: redshifted EPS; thin 
solid line: blueshifted EES. The observer is located at 2 — > 00. 



When integrating over the EFS one encounters a specific numerical problem related to 
sudden jumps of the brightness occurring at such pj(z/°°), at which the observer sees the 
boundary of the corresponding EFS(z/°°) (c.f. Figure EjlHl and also Paper I for discussion). 
If not treated properly, there are spurious oscillations of the line profile observed for certain 
sets of param eters. To some extent it ma y be cured by using finer discretization of the {pi} 
grid (see e.g. Marti fc NoerdlingeJ ( 1977 )). We prefer a different approach: for each z/°° we 
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Figure 3. Equal frequency surfaces; velocity law \24:) .(left). and l|25|l . with e = 0.1 {right). Other parameters: m = 2, go = 16, 
/3max = 0.15. Thick solid line: redshifted EPS; thin solid line: blueshifted EPS. 

numerically find such {zj,pj) at which dp/dz = 0, and then, spht the integral f l23|) . into two 
parts: 



which also allows to have two {pi} subgrids each having different mesh density. The latter is 
important for those cases in which different branches of EFS are tightly packed near compact 
object. Thus we insure that pj is always exactly on the boundary between cells and thus, 
we eliminate spurious oscillations of the line profile. 

6 LINE PROFILES 

For the given set of parameters, Qq, U°°, m, p, and e we calculate line profiles from equation 
(!23|) adopting velocity laws fj2^ . (!25|) . and the opacity law (126|) . All profiles are calculated 




(28) 
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assuming the size of the computational domain, Rout = 50 Re- In case of the failed wind law, 
(Elj), the value of initial velocity is fixed constant, wq = 0.01, and in all cases Vth = 57kms~^. 

6.1 Distorted P-Cygni and W-shaped profiles 

Figures m and Figure show profiles for different Qq, and velocity laws f l24p . and f l25p . From 
their corresponding upper left panels, we infer that both decelerating and failed winds show 
P-Cygni profile provided the influence of the gravitational redshifting is negligible. 

Notice that if equation ( 125|) is adopted, the parameter Uq does not equal the maximum 
velocity of the wind (see discussion after equation (l25l)). For example. Figures S] and [5] are 
calculated for U°° = Uq = 156.4, corresponding to V°° = 0.03c in the case of Figure SJ 
and Umax = 112, i.e. Wmax'^^th/c = 0.02 c in the case of Figure O At large 5^0 the maximum 
blueshift, i.e. y^ax — y°° — — 156.4, and the maximum projected velocity of the portion 
of the wind which is not obscured by the core determines the maximum observed redshifting: 
Vmax — ^67.6. Hereafter, is defined as a red- or blueshifted non-dimensional frequency 
y°°. For the decelerating wind Umax is reached very close to the bottom of the wind, where 
gravity is strong, opposite to the case considered in Paper I, where maximum velocity is 
reached far from the star where gravitational redshifting is negligible. 

The position of the red- ward edge of the line is determined from the approximate relation: 
?/iax - ^^max " (/(Sfi-o)- For = 25 wc havc = 52, and from upper right of Figure HJ 
we infer i/^ax — 46. Comparing to Figure H] we see that this approximate formula works well 
if "Umax is located close to the star. 

The velocity law flMj) generates the opacity distribution which has a maximum far from 
the star, while failed wind, ( 125|) leads to the opacity law, which peaks close to the pho- 
tosphere (c.f. Figured]). There are no W-shape profiles in Figure HI If opacity peaks close 
to the photosphere, as in Figure El then profiles with two absorption troughs separated by 
"emission" component are observed. 

From Figure [5] we can assume that W-shaped profiles are more likely to be observed from 
the failed wind (|25l) rather than from the case of the decelerated wind, ( l24l) . 

To verify this assumption we calculate profiles for the velocity law f l24l) adopting opacity 
distribution which peaks at the photosphere in case of the velocity law fl24|) : 

T,{w) = n{k + l)w\ (29) 

The results are shown in Figure |6l where both sets of profiles are calculated for qq = 10. 
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Figure 4. Line profiles for different launciiing radii. Optical depth law: I I26I I. with fc = 1; velocity law: I I24I I. with m = 1. Curves: 
dotted: To = 0.1, dashed: To = 0.4, dot-dashed: Tq = 1; solid: thin: Tq = 4; thick: To = 10. Other parameters: U°° = 156.4 
corresponding to V°° = 0.03 c. 



We conclude that a short acceleration phase of the failed wind is critical for the formation 
of the W-shaped profiles. 

One can notice small spikes of intensity apparent in Figure |6l they are most prominent 
in the absorption edge within the redshifted part of the profile (Figure [5], upper right). In 
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Figure 5. Line profiles from failed wind approximated by ||25|I . with m = 1 and e = 0.1; Uq = 156.4, which corresponds to 
"max-Uth/c = 0.02 c, ?^max — 112 (see text after 1I25I I about the difference between and Umax, and f/o)* Other parameters 
and the notation are the same as in Figure [4] 



Figure E] at the redshifted part of the spectrum they are observed for both of V°° = 0.03 c 
(left), and = 0.2c (right). The position of such a spike is independent of r and can 
be understood adopting the following considerations: If gravitational shifting is negligible, a 
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typical wind with f ~ 1/r law has EFS which looks as two symmetrical loops in p — z plane 
(for fixed ±?/°°) with respect of the z = plane (i.e. as in Figure [21 upper left). 

As Qo decreases, the redshifted loop covers the star from the observer, and at higher 
redshifts (lower go) the loop shifts further behind the star, and covers only a fraction of it 
(Figure El lower left). This is due to the combination of the gravitational redshifting and 
strong Doppler red-, blueshifting at the photosphere. For lower qq (higher redshifts) EFS is 
obscuring less of a stellar disk, and such partial coverage results in interception of a lesser 
fraction of the direct fiux of the core. Additional fiux in the spike comes from the continuum 
radiation scattered by the same EFS. 

The position of such spikes, y~ is approximately determined by gravitational and trans- 
verse Doppler redshift at the photosphere, y" ~ — C/(2(7o) — /3^C/2, at a; = 1; adopting 
parameters used to plot Figure |6] (left) we calculate y' ~ —262, while the figure gives 
y- = -291. 

The velocity at x ~ 1 is high and even a small deviation of /i from zero, produces 
significant normal Doppler blue- or redshifting; the effect is larger at larger v, e.g. for the 
case of V°° = 0.2 c (l6l right) we have: the predicted: ~ —364; the observed: y^ ^ —286, 
which is closer to the shift produced by the gravitational field alone, |/g^av — ~260. From 
the above considerations it is clear that position of such spikes is independent of the actual 
distribution of opacity. 

6.2 Inverted P-Cygni profiles 

Figure [7] shows profiles from failed and decelerated wind. In the case of the decelerated wind, 
varying m in ( l24l) or k in ( l26l) results in the distortion of the P-Cygni profile. If the wind 
decelerates more rapidly (larger m) then the system of EFS is located closer to the star, and 
the lower surface area of EFS results in lower intensity of the profile (upper left). The effect 
from increasing k in the opacity distribution is understood from Figure [11 i.e. in the case of 
the velocity law ( |2^ there is more opacity at the same radius, x if is larger, resulting in 
that both the local contribution, S{1 — e~^'), and the attenuation, /cC""^' is larger. 

For the same setup of parameters shapes of profiles from the failed wind are completely 
different. From Figure [7] (upper right) we see that as m increases, the P-Cygni profile first 
transforms to the W-shaped, and then to the inverted P-Cygni. 

To understand such behavior, we notice the W-shaped profiles are observed (at suffi- 
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Figure 6. Line profiles from decelerated wind with the velocity law I I24I1 . with m = 1, and opacity law I I29II : both sets of profiles 
are calculated for go = 10; = 0.03 c (left), V°° = 0.2 c (right). Curves: dotted; Tq = 0.1, dashed: Tq = 0.4, dot-dashed: 
To = 1; solid: thin: To = 4; thick: To = 10. 



ciently strong gravity) only in special cases of the accelerated wind (see Paper I), being very 
sensitive to the distribution of the opacity. 

In the case of accelerated wind in most cases there is a redshifted EFS behind the star at 
large x, because only Doppler redshifting is important there. In result the surface area of the 
corresponding redshifted EFS is larger and its emission overwhelms the redshifted absorp- 
tion. In the case of the decelerated or failed wind the maximum Doppler red-, blueshifting 
is made by the gas which is close to the star. 

To produce an absorption line the EFS should be in front of the star and needs to have 
surface area to cover enough of a star. Roughly the same shape, low surface area EFS is 
situated behind of the star, and i) if it has small surface area it reflects only a relatively 
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Figure 7. Line profiles for different velocity laws. Upper panels: curves for different m (i.e. slope of the velocity profile) in I l24t . 
and l|25ll . and fixed fc = 1 in l|26|l . Bottom panels: different k (i.e. power index in the opacity law) in lj26|l . and fixed m = 1 in 
| [24ll and | [25ll . Other parameters: go = 10, V°° = 0.1 c, Tq = 1, and e = 0.1; if the velocity law 1251 1 is adopted. 



small fraction of the flux towards the observer, and ii) because of the obscuration by the 
core the flux received by the observer is much lower. 

Thus, sufficiently strong gravity "quenches" the blueshifted EFS in front of the star, 
and thus there is no blueshifted absorption but only a strongly redshifted one. On both 
sides of the star Doppler blueshift from the gas is comparable with gravitational redshift 
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(e.g. roughly at fi < 7r/2). This emission adds to the radiation of the core. In result, the 
emission is blueshifted with respect to the redshifted absorption trough as shown in Figure [7] 
(upper right). Understanding that in the decelerated wind case the gravitational redshifting 
is capable of producing inverted P-Cygni profiles is one of the major findings of this work. 

It is instructive to investigate inverted P-Cygni profiles by changing the duration of the 
acceleration phase (i.e. e in (l25l)): the smaller e, the shorter is such phase (Figured]). In the 
case of the wind which is "almost" purely decelerating, e = 0.02 the results are shown in 
Figure ID (left). 

Increasing m has the following consequences: the distortion of P-Cygni profile (dashed); 
the formation of W-profile (dotted, thin solid), and finally the transformation to the inverted 
P-Cygni profile. The latter is characterized by a relatively weak absorption line superimposed 
on the blueshifted emission (thick solid). In the case of e = 0.2 the behavior is qualitatively 
the same: modified P-Cygni, W-shaped, and inverted P-Cygni. 

7 DISCUSSION 

Studying the influence of the gravitational redshifting on the shape of spectral line profiles 
we adopt two basic scenarios for the wind dynamics. In the Paper I we calculated profiles 
from monotonically accelerated winds (stellar-type winds, and homologous expansion). It 
was established that gravitational redshifting is of significant importance, being able to 
modify and distorting line profiles considerably. 

Profiles are sensitive not only to wind dynamics but also to the distribution of the opacity. 
Given the gravitational redshift is strong enough, and the terminal velocity, V°° is either 
much smaller than the escape velocity at the base of the wind, or the distribution of the 
opacity strongly peaks close to the star, then a new type of profile, a W-shaped profile is 
produced. 

We can imagine a situation (X-ray bursters as an example) when energy is deposited 
into the wind during a short initial pulse, and after that time the wind is decelerating. 
In this paper we consider winds which are either decelerating, according to v{r) ~ i/r™ 
law, or "failed", i.e. described by the bridging formula ( l25l) . In the latter case after a short 
accelerating phase, when v{r) ~ (r — -Re)™", the wind velocity decreases according to l/r™ 
law (c.f. Figured]). In both cases, the wind reaches its maximum velocity V°° in the region 
where gravity is strong. 
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Figure 8. Line profiles from failed wind lj29p : left: e = 0.02; right: e = 0.2. Curves are for different values of m at l)29p . Other 
parameters: go = 8i = 0.2 c. Opacity law II26I I with A; = 1. 



If there would be no absorption then the emission from the wind would be observed 
as a redshifted emission line having an approximately symmetric profile. In the case of the 
stellar-type wind, v{r) ~ (1 — Rc/r)°' far from the compact object Doppler blueshifting 
always overwhelms gravity, and as a result there is always a possibility to have blueshifted 
absorption line. This absorption eats away blue-ward part of the emission line and produces a 
P-Cygni-type profile. On the other hand, inner parts of the wind also have velocities which 
are not large enough for the Doppler blueshifting to overwhelm gravitational redshifting, 
and such absorption eats away the red-ward part of the emission line. The result is that in 
certain cases these reshifted and blueshifted lines are combined together being separated by 
the emission component producing a W-shape profile. 

The case of the decelerated wind is different. Such a wind has its maximum velocity from 
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the very beginning. The maximum Doppler blueshift is reached close to the photosphere, 
where sufficiently strong gravitational redshifting can overwhelm it. If there is an absorption 
it will be redshifted, and so it is basically the regime where the photospheric line eats away 
the red-ward part of the emission from the wind. In such a case the inverted P-Cygni profile is 
produced. Depending on how quickly the wind decelerates, there also may be an absorption 
occuring far form the star, i.e. eating away the blue-ward part of the line and producing a 
W-shape profile similar to the the case of the monotonically accelerated wind. 

Tree factors make inverted profiles possible: i)strong gravity which competes with the 
Doppler effect ii) the maximum velocity of the wind is approached in the region where 
gravitational redshifting is strong enough to overwhelm it, and iii) at larger radii the wind 
quickly decelerates. 

If the energy is deposited into the wind during a short initial pulse and the line is indeed 
formed in the vicinity of the compact object then it is likely to be transient, and profiles 
similar to those obtained in this paper should be considered as snapshots of different stages 
of expansion. The profile of a line from the early stage of expansion would be different from 
the profile of the same line at later times. Though we believe that obtained profiles are 
robust, the need of deducing physical conditions of the line forming plasma in a real object 
would question simple approximations adopted in this pa per, and such problem should b e 
treated using methods of the full 3D modeling (such as in 
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8 CONCLUSIONS 

Several conclusion can be drawn about the influence of strong gravitational redshifting on 
profiles of spectral lines from failed and decelerated winds: 

i) At large Rc/^g where gravitational redshifting is negligible, decelerated winds (!24l) and 
failed winds f l25|) produce profiles with shapes similar to each other and to P-Cygni proffie. 
Profiles from the failed wind have absorption troughs which are flatter than those from the 
decelerating wind. 

ii) If Rc is smaller than several tens of the effect of the gravitational shifting is significant. 
The strength of the effect depends on the duration of the acceleration phase (parameter e): 
profiles from decelerated wind (e = 0) are "equivalent" to P-Cygni, i.e. they have absorption 
trough which is blueshifted with respect to the emission peak. Changing the slope of the 
velocity proffie or changing the opacity distribution influences mostly the net intensity within 
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the emission/absorption components; if e > (failed wind) profiles are W-sliaped, sometimes 
with additional weak emission component at the red-ward part of the profile (saw-toothed 
profile). Such profiles are obtained in many cases when Rc is less than approximately 20- 
30 rg. The effect of the gravity does not depend on V°° as sensitively as in in the case of 
monotonically accelerated wind, 
iv) If failed wind is decelerating quickly (e.g. the parameter m > 3 — 4 in (!25l) . for = 
10 r^, V°° = 0.1c) the inverted P-Cygni profile is observed. Increasing the rate of the wind 
deceleration (i.e. m), first produces the distorted P-Cygni profile, then changes it to W- 
shaped, and then, at higher m produces the absorption-emission inverted P-Cygni profile. 
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